arXiv: 1504.01848vl [cond-mat.mes-hall] 8 Apr 2015 


Instability-induced formation and non-equilibrium dynamics of phase defects in 

polariton condensates. 

T. C. H. Liew, 1 O. A. Egorov, 2 M. Matuszewski, 3 O. Kyriienko, 1,4 X. Ma, 2 and E. A. Ostrovskaya 5 

1 Division of Physics and Applied Physics, Nanyang Technological University 637371, Singapore 
2 Institute of Condensed Matter Theory and Solid State Optics, Abbe Center of Photonics, 

Friedrich-Schiller-Universitdt Jena, Max-Wien-Platz 1, 077f 3 Jena, Germany 
3 Instytut Fizyki Polskiej Akademii Nauk, Aleja Lotnikow 32/f6, PL-02-668 Warsaw, Poland 
*QUANTOP, Danish Quantum Optics Center, Niels Bohr Institute, 

University of Copenhagen, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
5 Nonlinear Physics Centre, Research School of Physics and Engineering, 

The Australian National University, Canberra ACT 2601, Australia 
(Dated: January 23, 2015) 

We study, theoretically and numerically, the onset and development of modulational instability 
in an incoherently pumped spatially homogeneous polariton condensate. Within the framework of 
mean-field theory, we identify regimes of modulational instability in two cases: 1) Strong feedback 
between the condensate and reservoir, which may occur in scalar condensates, and 2) Parametric 
scattering in the presence of polarization splitting in spinor condensates. In both cases we investigate 
the instability induced textures in space and time including non-equilibrium dynamics of phase 
dislocations and vortices. In particular we discuss the mechanism of vortex destabilization and 
formation of spiraling waves. We also identify the presence of topological defects, which take the 
form of half-vortex pairs in the spinor case, giving an “eyelet” structure in intensity and dipole type 
structure in the spin polarization. In the modulationally stable parameter domains, we observe 
formation of the phase defects in the process of condensate formation from an initially spatially 
incoherent low-density state. In analogy to the Kibble-Zurek type scaling for nonequilibrium phase 
transitions, we find that the defect density scales with the pumping rate. 

PACS numbers: 71.36.+c,42.65.Sf,78.67.-n,03.75.Mn 


I. INTRODUCTION 

The formation of complex spatio-temporal patterns 
and textures is a particularly intriguing phenomenon oc¬ 
curring in a diverse range of physical systems!- 1 ^ The 
patterns commonly arise in nonlinear dissipative sys¬ 
tems driven far from equilibrium. Amongst examples 
of such systems, exciton-polaritons in semiconductor mi¬ 
crocavities have emerged as a hybrid light-matter sys¬ 
tem with strongly nonlinear properties. Effects such as 
Bose-Einstein condensation,®^ superfluidity®^ anc [ the 
formation of solitons^ 13 have been well-documented in 
the literature! 14 * 15 ! 

A particularly important nonlinear effect in the con¬ 
text of pattern formation in microcavities is the paramet¬ 
ric scattering of resonantly excited polaritons in planar 
semiconductor microcavities.®® The pair scattering of 
pairs of polaritons to different states in reciprocal space 
allows a homogeneous polariton field to spontaneously 
break translational symmetry. This effect enables the 
formation of ordered hexagonal/triangular lattices, pre¬ 
dicted theoretically®® an d observed experimentally— 
under different excitation conditions, as well as lattices 
of breathing solitons.® 

Under non-resonant/incoherent excitation, v ortex lat¬ 
tices were predicted to occur in harmonic trapd^® an d 
later observed in experiments involving multiple excita¬ 
tion spots!- 27 The formation of multi-lobec® and vortex- 
antivortex patterns 29 under ring shaped excitation, as 


well as sunflower rippled® excited by a narrow pump spot 
have also been reported. In these examples, the transla¬ 
tional symmetry of the system is already broken by the 
chosen shape of the pump-spot and/or the presence of a 
gradient in the potential of polaritons (as in, for example, 
the case of harmonic trapa^*). 

In this work we consider the possibility of spontaneous 
breaking of translational symmetry and pattern forma¬ 
tion in planar microcavities excited by a spatially homo¬ 
geneous incoherent pump. The pumping creates a reser¬ 
voir of “hot” exciton-like polaritons, which form a po¬ 
lariton condensate through a stimulated scattering pro¬ 
cess. The translational symmetry breaking is triggered 
by linear instability of the homogeneous condensate to 
spatial modulations, and the nonlinear evolution of the 
unstable state leads to formation of spatial patterns. We 
consider two different mechanisms of such modulational 
instability (MI) in this system. The first one arises when 
the polariton condensate has a strong feedback effect on 
the reservoir, in the form of reservoir depletion due to 
stimulated scattering of reservoir excitons into polari¬ 
tons. In this case, while polariton-polariton interactions 
are repulsive, the essentially saturable nature of exciton- 
polariton interactions may lead to effectively attractive 
nonlinearity for sufficiently low pump powers! 31 This ef¬ 
fective focusing nonlinearity in the system naturally leads 
to MI of a spatially homogeneous state, which was estab¬ 
lished in seve ral p revious studies for both quasi-ID and 
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Modulational instability is also known in spin-1 Bose- 
Einstein condensates of ultracold atoms due to paramet¬ 
ric couplin^2SH29| anc [ non ii near interactions between spin 
components,^ as well as in ID exciton-polariton conden¬ 
sates with a spin (polarization) degree of freedom.® Al¬ 
though polariton systems are non-conservative and non- 
equilibrium, the two-component spin degree of freedom 
of polaritons does allow a second mechanism of MI, 
which works also in the defocusing regime where a strong 
condensate-reservoir feedback is unnecessary. A circu¬ 
larly polarized excitation splits the energ y of th e cr + and 
o — states due to anisotropic interaction d 42 ! 43 ! occurring 
between the spin polarized reservoir and condensed po¬ 
laritons. This splitting sets the foundation for a paramet¬ 
ric scattering process as polaritons in an initially homoge¬ 
neous state with wavevector k = 0 on the upper spin-split 
branch can now scatter to degenerate non-zero wavevec¬ 
tor states on the lower branch, reminiscent of experi¬ 
ments under resonant excitation in triple microcavitied® 
or experiments in one-dimensional polariton systems.® 
While such a scattering process is not strictly allowed in 
isotropic cavities as it would violate spin conservation, 
the presence of sample anisotropy, which typically causes 
an additional linear polarization splitting and hybridiza¬ 
tion of the cr + and <r_ branches, relaxes this limitation. 

By considering the stability of the steady states of the 
system to weak perturbations, we find the zones of MI 
in the two different regimes. In the scalar case, where 
MI is derived from the condensate-reservoir feedback, we 
find that the homogeneous state breaks its translational 
symmetry and forms a turbulent state of phase disloca¬ 
tions, i.e., vortices. Unlike the previously studied cases, 
vortices do not form as the result of thermal fluctua- 
tiond® or scattering on disorder®. Rather, the spatial 
fragmentation of the initially homogeneous condensate 
due to the development of MI creates multiple interfer¬ 
ence between polariton flows generated by the randomly 
distributed sources, which leads to the development of 
multiple phase dislocations, similar to the scenario pre¬ 
viously considered for multiple pump spotd® and highly 
inhomogeneous trapped polariton condensate d 26 * 49 ! 

In the case of modulationally stable background, we 
show that multiple phase singularities can appear as a 
result of mean-field evolution of an initial white noise 
state, which mimics a pre-condensate state lacking spa¬ 
tial and phase coherence. Remarkably, formation of mul¬ 
tiple vortices in this scenario seems to be analogous to, 
but not the same as, the Kibble-Zurek mechanism, which 
acts during the quench through a phase transition to the 
Bose-Einstein condensation (BEC(p2ED. Indeed, the later 
describes the formation of boundary defects between dif¬ 
ferent domains of condensate which develop an indepen¬ 
dent phase rather t han inheriting it from the neighbour¬ 
ing spatial domaind- 52 * 53 ! 

We stress that the process of defect formation dur¬ 
ing nonequilibrium condensation of exciton-polaritons 
does not fo llow the scenario of the Kibble-Zurek mecha- 
nisn i 51 l 5 4 ^. The main difference is that in the latter, it 


is assumed that the system is initially in thermal equilib¬ 
rium, and is driven out of equil ibrium only in the vicin¬ 
ity of the phase transitior l 51 l 55 [ The process is divided 
into three phases, corresponding to adiabatic-impulse- 
adiabatic evolution. In nonequilibrium condensation, the 
system is far from equilibrium at the outset, and the tran¬ 
sition to the quasi-equilibrium (condensed) state occurs 
only after crossing the critical point. 

Nevertheless, the Kibble-Zurek mechanism and the de¬ 
fect formation in nonequilibrium systems have much in 
common. In both cases, defects are created due to sym¬ 
metry breaking in separate parts of the system which 
cannot communicate in a finite time. In both cases, there 
is a competition of two timescales existing in the system, 
which results in the same algebraic forms of power-law 
scalings for the number of defects and their characteristic 
creation timr®. In the polariton condensation case, the 
quench time is replaced by the timescale of the forma¬ 
tion of the condensate, which is controlled by the exter¬ 
nal pumping rate. We refer the reader to Sec. |IVD| and 
Ref. [HU for the detailed description of the process. 

Regardless of the mechanism of the vortex formation, 
either as the result of the MI development or as a result 
of transition to BEC, we show that the presence of the 
incoherent reservoir affects substantially both stability of 
vortices and their collective dynamics even for the case 
of a stable homogeneous background. As a consequence, 
the vortices can lose their radial symmetry and develop 
either into spatially localized rotating phase dislocations 
or into non-localized spiraling waves. 

A similar situation occurs in the spinor case, although 
multiple branches of modulationally unstable and stable 
solutions are present. Defects in the spin polarization 
of the condensate may appear even in the modulation- 
aly stable regime. Such structures move randomly in 
the microcavity plane and are composed of half-vortex® 
half-antivortex pairs,® exhibiting an associated dipole 
type spin texture. We predict that the density of vor¬ 
tices grows with increasing pump power similarly to the 
Kibble-Zurek scaling behaviour. 

The paper is organised as follows. In Sec. [TTJ we de¬ 
scribe the mathematical model of a semiconductor mi¬ 
crocavity operating in the strong coupling regime un¬ 
der incoherent homogeneous optical pump of a circular 
polarization. Then, in Sec. m we study the stability 
and collective dynamics of phase dislocations in a single¬ 
component polariton condensate. Here the dynamics is 
mostly affected by the modulational instability originat¬ 
ing from the strong feedback between the condensate and 
reservoir. In Sec. m we report a numerical analysis of 
the condensate dynamics in the presence of polarization 
splitting in spinor condensates. In Sec. |IVD[ we study 
the defect formation and the scaling laws for their density 
in analogy to the Kibble-Zurek mechanism. 
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II. THEORETICAL MODEL 


Let us begin by considering the incoherent excitation of a 
spinor polariton system, i.e., the system where the polar¬ 
ization degree of freedom is significant. The scalar case, 
which is valid when only one spin component is pop¬ 
ulated, is then easily obtained by removing one of the 
spin components. It is worth recalling that experiments 
with a circularly polarized optical pump have resulted 
in the excitation of a circularly polarized polariton con¬ 
densate at the pump position, in both 2D^ and ID— 
samples. A theoretical model can be based on the gener¬ 
alized Gross-Pitaevskii approach,^ where a condensate 
of exciton-polaritons can be described by the wavefunc- 
tions, ip+ and ip-, of the a + and cr” circularly polarized 
states, respectively: 

= (~”2uT +9lTlR + a il^+! 2 + a M-\ 2 ^j ^+ 

+ ih (rn R - T) ip + + Axy^ 1 -, ( 1 ) 


. t dip- ( H 2 \/ 2 . . l2 . . l2 \ . 

= l-+ g 2 n R + a 2 \tp+\-+ axl'ip-l \ ip- 

— ihTip- + AxyVM-• (2) 

Here we assume that the circularly polarized pumping 
creates a circularly polarized reservoir, n R , with dynam¬ 
ics described by the rate equation: 

= - (r fl . + r\ip + \ 2 ) n R + P, (3) 

where P represents the pumping rate. Non-linear 
interactions between polaritons are characterized by a± 
and 02 , representing the int eracti on strengths between 
parallel and antiparallel spin d^ 2 ^ , respectively. Similar 
parameters, g\ and 92 , characterize the blueshift caused 
by the circularly polarized reservoir. T is the polariton 
decay rate. A y y represents a linear polarization split¬ 
ting, which has been reported in several ex perim ental 
studiedMIiH and can take values of 0-0.2 meV “2122] Even 
larger values can be expected by the application of 
magnetic fields (in Voight configuration). 


We note that in the limit T# T, it is possi¬ 
ble to pr oceed by adiabatic elimination of the reservoir 
dynamics! 26 * 64 ! However we do not make such an approx¬ 
imation here since the reservoir dynamics is important 
for our further analysis. 


III. NON-EQUILIBRIUM DYNAMICS IN THE 
SCALAR CASE 

A. Modulational instability of the homogeneous 
steady state. 


If we consider a circularly polarized pump in a micro¬ 
cavity with negligible polarization splitting (Axy = 0) 




FIG. 1: (Color online) Steady-state solutions of the model (JTJ- 
and their stability for a circularly polarized pump (Axy = 
0). (a) The homogeneous solution (HS) (|V , o| 2 ) as functions 
of pump P. The dotted line depicts the modulational in¬ 
stability (MI), (b) Growth rates [Amu; (ps -1 )] of the small 
perturbations (around |V'o| 2 ) as a function of momentum k 
and HS density |'i/’o| 2 - The vertical axis is the same as in (a) 
and the horizontal dashed (blue) line connecting (a) and (b) 
indicates that the boundary of modulational instability oc¬ 
curs for the same maximum homogeneous density. Snapshots 
of the condensate phase in real space (c) and the density in 
two-dimensional reciprocal space (d) within the MI domain of 
a HS for pump P = 9 ps ~ 1 gm~ 2 . k cu t 0 a can be identified as 
the magnitude of the in-plane wavevector at which all unsta¬ 
ble modes disappear. Parameters: ai = 6 x 10 -4 meVgm 2 , 
gi = 4ai, T = 0.165 ps -1 , r = 0.01 ps -1 /xm 2 , Fr = 0.5 ps” 1 . 
The polariton effective mass was taken as x 10 -4 of the free 
electron mass. 


then the model can be reduced to a scalar one (ip = ip+) 
without the second polarization component (ip- = 0). 
First, we study a scalar steady-state homogeneous solu¬ 
tion (HS) of the system of Eqs. 0-© and discuss its 
stability. For the sake of generality we allow also the 
HSs with nonzero transversal momenta k 0 ^ 0, which 
have the form of travelling waves 

ip hs (x,t) = ipoe-^V* ( 4 ) 

where the condensate energy is given by h/j, (ko,\ipo\) = 
(h 2 /2m) k,Q + ai\ipo\ 2 + giriRo. The HS becomes nontriv¬ 
ial provided that the external pump compensates for all 
losses and overcomes the threshold value P^l P th = TT R /r. 
The coherent exciton-polariton density and incoherent 
reservoir density are given by \ipo\ 2 = (P — Pt/ 1 )/T and 
n R0 = r/r, respectively. 

The linear stability analysis of the homogeneous steady 
state of our scalar system and its modifications has been 
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previously performed by many authors! 15 ! 31 * ^ For our 
choice of the system parameters, the linear stability anal¬ 
ysis shows that the HS becomes modulationally (dynam¬ 
ically) unstable within a pump interval just above the 
threshold value of the pump P t h [Fig. [lja)] (details of 
the analysis are given in Appendix [A}. This MI is asso¬ 
ciated with the parametrical generation of field compo¬ 
nents with nonzero momenta k. Figure 0 b) presents the 
linear growth rate Smw(fc) > 0 of the unstable pertur¬ 
bations as a function of their momenta k. 

It has been shown recently®] that this MI is associ¬ 
ated with the effective attractive nonlinearity induced 
by the saturation of the incoherent reservoir. Indeed, 
based on the intuition gained from paradigm nonlinear 
models, such as the Schrodinger equation with a Kerr- 
nonlinearity, one expects that the existence of the MI 
requires a focusing nonlinearity.®^ 68 However, owing to 
repulsive interactions between excitons, the nonlinear be¬ 
haviour of an exciton-polariton condensate is akin to 
that of optical waves in a defocusing media. This seem¬ 
ing contradiction clearly elucidates the influence of the 
open-dissipative nature of the system on the nonlinear 
behaviour,® which requires inclusion of an incoherent 
reservoir of “hot” excitons. To illustrate this influence 
we consider a nonlinear energy shift induced by both the 
coherent exciton-polaritons and the incoherent reservoir 


hum (|V’| 2 ) = <*iM 2 +5i-- P , ,, 2 - ( 5 ) 

v ' T R + r\ip\ 

We note that the polariton density \ip\ 2 corresponds to a 
steady-state solution which, in general, is not necessar¬ 
ily given by the homogeneous value IV’ol 2 - The reservoir 
intensity n R follows this steady-state solution \ip \ 2 . The 
first term of Eq. ([5]) describes the blue shift originating 
from the repulsive exciton-exciton interaction, whereas 
the second term describes the saturation of the reser¬ 
voir. We define the effective nonlinearity coefficient as 
g e g = hdf-i n i/d\ip \ 2 . Then the nonlinear response is ef¬ 
fectively focusing provided that the coefficient is nega¬ 
tive g e g < 0. Otherwise, if g e g > 0 the nonlinear re¬ 
sponse is defocusing or repulsive [see thick solid lines in 
the Fig. [2](a)] . In the vicinity of the steady-state HS with 
intensity |i/’o| 2 this coefficient takes the form (cf. Eq. (36) 
in Ref.[SI]): 

«=»■ 0 - S 5 ) ■ (e) 

The effective nonlinear coefficient [Eq. ([6])] changes sign 
for a pump value Pmi = giT 2 /air. It means that in the 
vicinity of the HS the nonlinear response changes charac¬ 
ter from effectively focusing to effectively defocusing. As 
a result the modulationally unstable HS becomes stable 
exactly at this point P = Pmi- The condition Pmi > Pth 
gives a general criterion for the appearance of MI derived 
in Ref. [2T1 : 


9 iT 


> 1 . 


(7) 




FIG. 2: (Color online) (a) Effective nonlinear energy shift 
hg n i versus the condensate density, \if>\ 2 , for different pump P. 
Dashed lines (with the negative slope) represent an effectively 
focusing nonlinear response (or nonlinear red shift). The grey 
triangle represents the intensities \ip\ 2 exceeding values given 
by the steady-state homogeneous solutions (|i/’o| 2 )- (b) Vortex 
width (FWHM) vs. the pump P. The dashed line represents 
the width approximated by the value ss 3.5£. d max depicts the 
maximal diameter allowed for stable vortices, see Eq. (<0 
The radial dependence of the condensate and reservoir density 
within the vortex profile for P = 20 ps -1 pm~ 2 . (d) The real 
and imaginary parts of the tp within the vortex profile given 
in (c). Other parameters are the same as in Fig. [T] 


Within the MI interval ( P t h < P < Pmi) indicated in 
Fig. 0a), the HS experiences spontaneous translational 
symmetry breaking, resulting in the formation of non- 
uniform turbulent states of the condensate. 

The result of a direct numerical calculation is shown 
in Fig. 0c). In these calculations and those presented 
throughout the manuscript we use a square grid with at 
least 128 x 128 points covering the plot area. We make 
use of an adaptive step Adams-Bashforth-Moulton pro¬ 
cedure, which was previously found to be consistent with 
fixed step methods with an integration step of O.Olps. 
To double check the results of our numerical simula¬ 
tions we repeated them with higher precision using up 
to 1024 x 1024 grid points. Unless stated otherwise, pe¬ 
riodic boundary conditions are applied. 

The strongly non-equilibrium state in Fig. 0c) in¬ 
cludes one- and two-dimensional phase dislocations which 
move chaotically and overlap. Therefore this dynamics 
can be characterized as a “strong” turbulence regime 
with overlapping defects (see Ref. [48] and references 
therein). Two-dimensional Fourier transformation shows 
that the most part of the spatial momenta are bounded 
within the ring with the radius given by the cutoff mo¬ 
menta of the unstable modes fc cuto ff [Fig. [l|(d) l . N ote 
that direct numerical simulations of the model (|1|3[) for 
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different initial conditions did not reveal the formation of 
stationary periodical patterns, known for the coherently 
pumped polaritonic systems®®] This is due to desta¬ 
bilization caused by fluctuations in the exciton reservoir. 
The modulational instability discussed later in the spinor 
case (section IV) does not depend on having a dynamic 
reservoir and can largely be reproduced assuming a static 
reservoir2S. This allows an explicit testing of the effect of 
the dynamic reservoir in the spinor case, where we find 
that periodic patterns are possible with a static reser¬ 
voir but are prevented as soon as the reservoir density is 
allowed to evolve spatiotemporally. We assume that it 
is the freedom for density fluctuations to appear in the 
reservoir that lead to the disruption of regular patterns 
in the condensate also in the scalar case. 


B. Single vortices in the dynamically stable regime. 


the effective focusing nonlinearity (fl'effGV'l ) < 0) is al¬ 
ways satisfied within the vortex profile close enough to its 
core. Therefore the steady-state condensate on a circle 
with a fixed radius g c = \Jx 1 + y 2 and density \ip{g c )\ 2 
can become unstable against spatially modulated per¬ 
turbations. The periodical boundary conditions on the 
circle restrict the number of available momenta to the 
values k c — n/g c , where n is an integer. In the limit of 
large radius g c —> oo the values of the unstable momenta 
can be approximated by those values calculated for the 
HSs [see Fig. [ljb)]. Therefore the polariton condensate 
becomes unstable provided that at least one of the avail¬ 
able momenta ( k c ~ n/g c ) is smaller than the cut-off 
value fc cu toff ~ 0.19/rm -1 for MI. This allows estimation 
of the vortex diameter d max where the instability just 
sets in, i.e., fc cu toff = 2 /d max . Therefore the vortices are 
stable provided that their diameters do not exceed the 
maximal value given by 


The reservoir contribution becomes negligible in the 
limit of very strong pump P P t h■ We expect that in 
this case the nonlinear dynamics is very similar to that 
known for the conservative systems, including the forma¬ 
tion of the stable phase dislocations and vortices. Indeed 
the scalar version of the equations possesses vortex solu¬ 
tions within the stability interval of the HS for P > Pmi 
[see Figs.[2](b) and (c)]. Similar to the conservative case, 
one can define a characteristic length or an effective heal¬ 
ing length in the vicinity of the HS \ifo\ 2 '- 


dmax — 2fc cuto ff- (^) 

It is shown in Fig. [2](b) that the vortex diameter exceeds 
the maximal value d max within the pump interval Pmi < 
P < 18 ps^Vm" 2 in the vicinity of the MI do main. 
As we will see below (in the next subsection III CI they 
become unstable and lose their radial symmetry. 


C. Collective dynamics of vortices 


m = , ( 8 ) 

y2m <7eff|V'o(-P)| 2 

The healing length is a typical length scale over which if 
can change significantly. It also gives the typical size of 
the vortices. More precisely the full width at half maxi¬ 
mum (FWHM) of the vortices is approximately given by 
d « 3.5£ at least for a strong enough pumping [P Pth 
in Fig.|];b)]. It is worth to mention that the spatial oscil¬ 
lations of the real and imaginary parts of the amplitude 
profiles of the vortices go substantially beyond the heal¬ 
ing length [Fig. |2][d)], especially for a moderate pumping 
in the vicinity of the MI. This indicates permanent en¬ 
ergy and polariton exchange within the vortex profile. 
The presence of these intrinsic fluxes essentially changes 
the interaction dynamics between vortices!^ providing a 
purely dissipative mechanism for their mutual repulsion. 

The radial symmetry of the vortex phase is broken also 
for a n inh omogeneous pump, for instance, a Gaussian 
pumjjHEID As a result, the vortices have spiraling phases 
indicating again the permanent exchange of particles be¬ 
tween different points within the resonator plane. 

The stability analysis, discussed in the context of the 
HSs, can also be applied to the vortices, at least in the 
limit of very broad states. Indeed, following the heal¬ 
ing length £(P), the vortex width increases for smaller 
values of the pump and diverges in the vicinity of the 
MI [P = Pmi in Fig. [2jb)] . The formal condition for 


Beyond the MI instability interval for P > Pmi the 
non-trivial HS solution is stable. However, at the onset of 
condensation when it passes into the mean-field regime, 
the polariton state is spatially incoherent, with random 
phase. This leads to the generation of phase defects as 
the system goes through the condensate transition, in 
analogy to the Kibble-Zurek theory. To go beyond the 
mean-held approximation and describe first and second 
order spatial coherences when crossing the condensation 
threshold, one can make use of stochastic classical held 
approachesP'OEU These describe the pre-condensate state 
as an ensemble of fluctuating white noise states (governed 
by a stochastic Gross-Pitaevskii equation). The projec¬ 
tion onto the classical condensate state upon condensa¬ 
tion selects a particular realization of the noise. After 
the condensate has formed we assume that huctuations 
are weak in comparison to the condensate mean-held and 
can be neglected. This approach reproduces the typical 
establishment and coherent evolution of topol ogical de¬ 
fects in polariton condensate experiments®^!. We note 
that the account of huctuations throughout the evolution 
would be important for describing accurately spectral po¬ 
lariton properties or the polariton photoluminescence be¬ 
low threshold' 0 . Fluctuations can in principle shift the 
phase boundaries between stable and unstable region^?-, 
however, these shifts are expected to be limited and have 
little further effect on the dynamics. 

For a pumping substantially above the MI threshold 
the dynamics is dominated by the defocusing nonlinear- 
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FIG. 3: Nonlinear dynamics of the condensate beyond the 
MI threshold for (a,b,c) P = 25 ps -1 /r m~ 2 and (d,e,f) 
P = 17.5 ps ~ 1 pm~ 2 . (a,d) Long-time evolution dynamics of 
the dislocations number. (b,e) Snapshots of the condensate 
density (pm~ 2 ) profiles. (c,f) Snapshots of the condensate 
phase profiles. See also the Supplemental Material^ 74 ! Pa¬ 
rameters are the same as in Fig. [I] 


ity of a pure condensate (a). The influence of dissipa¬ 
tive dynamics of the reservoir becomes less important. 
Using the scalar version of the governing equations ([l]), 
© we calculated the condensate dynamics starting with 
a spatially incoherent state given by a small-amplitude 
white noise. Our initial condition mimics a particular 
realisation of the stochastic pre-condensate state. Sim¬ 
ilar to the conservative limit we observed the forma¬ 
tion of a spatially coherent condensate accompanied by 
the spontaneous formation of the vortices which move 
chaotically and interact with each other. Two vortices 
with equal (opposite) topological charges repel (attract) 
each other. Two attracting vortices mutually annihi¬ 
late if the distance between them becomes smaller than 
the healing length. Therefore, as discussed in Ref. [48] , 
the whole number of vortices drops gradually with time 
and approaches zero, at least for a very strong pump 
P > 40 ps ~ 1 pm~ 2 . 

In contrast, for a weaker pump, the number of phase 
singularities converges eventually to some constant value 
indicating the formation of a coherent state with a fi¬ 
nite number of dislocations, i.e., a superfluid turbu¬ 
lence [Fig. [3](a)]. A snapshot of the intensity and phase 
profiles shows a state of well distinguishable vortices 
[Figs. ©b) ,(c)]. The vortices move chaotically, inter¬ 
acting with their neighbours, and, in general, sustain a 


dynamical equilibrium. It is remarkable that the aver¬ 
age separations between nearest vortices remains more 
or less constant for this particular interval of pump val¬ 
ues. This means that there exists some equilibrium dis¬ 
tance between vortices. Apparently the influence of the 
dissipative effects and the conde nsate flows (mentioned 
in the previous subsection III B ) are substantial.^ The 
out-going condensate flows from the vortex centers hinder 
attraction between vortices and their annihilation. For 
even weaker pump the dissipative effects become stronger 
and, as a consequence, the average distance between vor¬ 
tices under dynamical equilibrium becomes even smaller 
[Figs. [3](d),(e) and (f)]. We note that this dynamical 
equilibrium forms over a long time scale exceeding hun¬ 
dreds of nanoseconds. Therefore the Kibble-Zurek-like 
scaling law does not describe the number of vortices in 
this regime (see Sec. IVD below). 

Even though the HSs are stable, the nonlinear dynam¬ 
ics of vortices is strongly affected by the reservoir sat¬ 
uration dynamics. The numerical modelling shows that 
the vortices themselves become unstable and develop into 
radially asymmetric rotating structures, as can be seen 
on the snapshot profiles in Figs. [3])e),(f) (a movie show¬ 
ing the time dynamics is available in the Supplemental 
Material 74 ]). This is in agreement with the destabiliza¬ 
tion scenario for the vortices discussed in the previous 
subsection. Indeed the vortex size exceeds the maximal 
diameter d max given by Eq. © and therefore becomes 
unstable. 


D. Formation of spiraling waves 

In the vicinity of the MI the vortices are strongly unsta¬ 
ble and the initial noise first develops into non-uniform 
dynamical states similar to those which appear for the 
modulationally unstable background [Figs. S a >b)]. In 
this “strong” turbulence regime® the characteristic dis¬ 
tance between vortices is substantially smaller than their 
typical core size, so the vortices are not structured and 
the chaotic behavior is seen on the level of a single 
vortex. However, after some sufficiently long time of 
about several thousand polariton lifetimes, the system 
switches spontaneously into a more regular regime, char¬ 
acterized by the formation of a single spiraling topo¬ 
logical dislocation [Fig. [4] (a movie showing the time 
dynamics is available in the Supplemental Material^). 
This spiraling topological state drives away other phase 
dislocations in the system and covers eventually the 
whole computational window provided that Pm t < P < 
16 ps ~ l pm~ 2 . Similar spiraling waves are known for 
other non-equilibrium dissipative systems! 75 ' 6 In the 
simplest case they are solutions of the complex Ginzburg- 
Landau equation. 

It is worth to articulate the differences between the spi¬ 
raling topological states and the vortices. First, there is 
no rotational symmetry in the profiles [Figs. 4]. Second, 
these topological states are not stationary and experience 
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top” shape in the form of a super-Gaussian intensity dis¬ 
tribution. It turned out that the spiraling waves appear 
also for the localized pump where the condensate density 
vanishes at the boundaries of the computational window. 
These calculations serve as a solid proof of the existence 
of the spiraling waves independently of a particular choice 
of the numerical boundary conditions. 

In general the out-going radiation from the center of 
the topological solution repels the local inhomogeneities 
of the profile and other topological solutions. This gives 
an additional purely dissipative mechanism which en¬ 
forces a long range spatial coherence in non-equilibrium 
systems operating in the regime of “strong” turbulence. 


IV. NON-EQUILIBRIUM DYNAMICS IN THE 
SPINOR CASE 


FIG. 4: (Color online) Snapshots of the condensate intensi¬ 
ties (a,c,e) and phase profiles (b,d,f) for different time points 
and P = 13.5 ps"^mr 2 . (a,b) Initial turbulent state of the 
condensate at t = 1000 ps. (c,d) Onset of a spiraling wave at 
t = 5000 ps. (e,f) The spiraling waves at t = llOOOps. The 
rotation period of the spiraling wave is 144ps. See also the 
Supplemental Materials' 4 . Other parameters are the same as 
in Fig. [T] 


uniform rotation of the density with a typical rotation pe¬ 
riod of about 144ps. Third, far from the center the pro¬ 
file converges to the homogeneous traveling wave solution 
given by Eq. Q with the amplitude ipo and a nontrivial 
momentum k o. Apparently this solution is characterized 
by a permanent radial flux of exciton-polaritons from the 
vortex center to the periphery and therefore resembles a 
point-like source of ring waves. We note that these spi¬ 
raling waves can exist only in nonequilibrium dissipative 
systems. 

An important question remains whether the spiraling 
state with a non-zero orbital angular momentum can 
emerge from the initially non-rotating turbulent state. 
Indeed, according to the conservation law of the total or¬ 
bital angular momentum, the phase dislocations appear 
in pairs which is also valid for the turbulent state con¬ 
sidered here. Apparently the local intensity fluctuations 
break the symmetry between the two dislocations within 
a pair in such a way that only one of them develops into 
the spiraling wave. The periodical boundary conditions 
(in x and y directions) were used in our numerical mod¬ 
elling. However, we have confirmed that the spiraling 
waves appear with the same probability independently 
on the computational window size and particular reali¬ 
sations of the initial seeding noise. Moreover, since some 
phase dislocations are always present [see Figs. |4je) and 
(f)], the total orbital angular momentum of the conden¬ 
sate within the computational window remains zero. We 
performed additional numerical simulations of the con¬ 
densate dynamics under a localized pump with a “flat¬ 


A. Stability of Homogeneous States 


In the presence of non-zero polarization splitting 
(Axy 7^ 0), the population of a~ component appears 
even for a fully polarized a + pumping. To emphasize 
the difference to the nonlinear dynamics discussed above 


(see Sec. Ill) we consider system parameters which do 
not satisfy the MI criteria ([7| and, therefore, guarantee 
the stability of homogeneous solutions in the scalar limit. 

Homogeneous stationary solutions can be found by 
substituting trial solutions in the form ip± (t) = ip±e~ lllt 
into Eqs. 0 and ([2]). This gives four stationary equa¬ 
tions for real and imaginary parts of the amplitudes. 
These are supplemented by the requirement that the time 
derivative in Eq. ([3]) vanishes for a steady state, that is, 

\ip+\ 2 = y — T^j. Noting that the phase reference 


of the system can be freely chosen, we can for simplicity 
set = 0, which allows one to find the relation 

$ie{ip_} = A p v 8m from one of the four station¬ 

ary equations. The remaining three equations can be 
solved for the remaining unknown quantities: Jfe {'!/’+}, 
3m {'0+} and p. 

The dependence of the stationary HSs on the pumping 
power is shown in Figs. 5ja) and (b). Linear stability 
of the steady states can be determined by the standard 
extension of the Bogoliubov-de Gennes analysi^^A ti — ^ 
by considering perturbations in the polariton and reser¬ 
voir fields of the form 5^ ± = U±e i(kx-ut) +v *± e -i(kx-iuj*t) 

and Sn R = w , respectively. The 

details of the derivation can be found in Appendi x [A} 

The stationary states labelled U\ and U 2 in Fig. f5ja,b) 
are unstable to fluctuations with k = 0. This “single¬ 
mode instability” indicates that even in a confined system 
(e.g., micropillar) the homogeneous steady state would be 
linearly unstable in the Lyapunov sense and any spatially 
homogeneous perturbation would grow. The stationary 
states labelled MIi and ML are unstable to spatial mod¬ 
ulations with non-zero wavevectors k. These states would 
be stable in a confined system, with no spatial degrees 
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FIG. 5: (Color online) Stationary solutions to Eqs. m w 
and (b) show the dependence of the a+ and <r_ polarized 
polariton populations on the pumping strength, respectively. 
The different branches are labelled stable (S), MI, and single¬ 
mode unstable (U). (c) Parametrically unstable zones of the 
MIi and MI 2 branches. The shading illustrates the gain of 
the instability (given by the imaginary part of to) for the MI 2 
solution. The grey rectangle denotes a region below conden¬ 
sation threshold where 1 p± = 0. The vertical line represents 
the pump intensity considered in Figs. [7] and Fig. [9] Param¬ 
eters: ai = 1.55 x 10~ 4 meV/im 2 , 0:2 = — O.lai, gi = an, 
(?2 = « 2 , r = 0.033 ps-\ r = 0.01 ps _1 /rm 2 , Tu = 10 ps _1 , 
A xy = 0.1 meV, P = 92 ps _1 pnG 2 . The polariton effective 
mass was taken as 5 x 10~° of the free electron mass. 


of freedom, however, in a 2D planar system the states 
MIi and MI 2 undergo parametric scattering. The state 
labelled S is fully stable, as the imaginary part of E = frio 
remains negative for all wavevectors (see Appendix 0 - 

The “S”-shape of the curves MI 1; MI 2 and U 2 is char¬ 
acteristic of multistability, which is a common feature 
of resonantly excited microcavitie^Ksni but less studied 
under the non-resonant or incoherent excitation 8 ^ that 
we consider here. While multistability is strictly only 
present in the confined system, since the MIi and MI 2 
states are unstable in the presence of spatial degrees 
of freedom, they can still give rise to different (non¬ 
stationary) configurations under the same excitation con¬ 
ditions. 

Fig. E c) shows the regions of MI in the system when 
the MIi or MI 2 branches are excited. The MIi branch 
begins once the threshold for polariton condensation is 
passed (indicated by the grey rectangle), such that only 
weak pump intensities are needed to see the effects of MI. 


B. Spin Textures due to Modulational Instability 

Due to the presence of MI, we can expect the frag¬ 
mentation of the homogeneous density of the condensate 
and spontaneous formation of spin textures, even in the 
presence of homogeneous pumping. Solving Eqs. 0-© 
numerically when the system is excited just above thresh¬ 
old on the MIi branch, we obtain the spin texture shown 
in Fig. [6] In analogy to the scalar case, the texture com¬ 
prises a simultaneous modulation of the intensity and 
phase in the system, which oscillate with multiple fre¬ 


quencies. In addition to the scalar case, there is also 
an appearance of a non-uniform polarization, despite the 
fact that the pumping of the system is homogeneous in 
both intensity and polarization. 
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FIG. 6: (Color online) Spinor MI regime, (a) Density pat¬ 
tern of a polariton condensate under homogeneous incoherent 
pumping, (b) Distribution of the circular polarization degree, 
pz = (IV'+I 2 - IV'-I 2 ) / (IV'+I 2 + l'0-l 2 )- (c) phase of the ip+ 
polariton field (the i/)_ component (not shown) has a similar 
dependence), (d) Time evolution of the polariton intensity 
along a slice in real space. Parameters were the same as in 
Fig. 0 (a small pump power was chosen so as to excite only 
the MIi branch). 


C. Spin Defects in the Dynamically Stable Regime 


When the system is excited with a larger pump power, 
the system tends to follow the stable branch (S) in Fig. [ 5 ] 
In this case one can expect a spatially HS due to the sta¬ 
bility of the (S) branch, however, defects present in the 
initial state after transition to condensation (in simula¬ 
tions taken as a low intensity white noise as in sections 
III C and III D) are trapped in the system and stabilize 
with the structure shown in Fig. [7] 

The defects have a non-trivial density structure with 
localized maxima inside of an otherwise circular shaped 
drop in density. These “eyelets” are not fixed in their 
locations, but move randomly in the plane with a typi¬ 
cal speed on the order of 0.05 /xm/ps. The eyelets also 
possess a characteristic spin polarization with a dipole 
type shape, as shown in Fig. 7](b). The structure of the 
eyelets can also be understood when looking at the phase 
distribution, which is shown in Figs. [Tfc) and (d) for the 
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FIG. 7: (Color online) Spin defect formation without MI. a) 
Density pattern of a polariton condensate under homogeneous 
incoherent pumping, b) Distribution of the circular polar¬ 
ization degree, p z = (\ip+\ 2 - |V'-| 2 ) / (IV , + | 2 + IV , -| 2 )- The 
black curves show the path traced out by the eyelets over 
1000 ps. c) and d) show the phase of the ip+ and if)- polari¬ 
ton fields, respectively. Black spots indicate the positions of 
(half) vortices. Parameters were the same as in Fig. [5] (with 
P = 84 ps - Vm -2 ). 


cr + and er _ polarized components, respectively. Here it 
is clear that each eyelet is formed from a pair of half 

vortice^Pl - 

one appearing in each spin component. 

In addition to the slow drift, the eyelets preserve their 
shape while undergoing a faster periodic rotation with a 
period close to 2tt/Axy (a movie of the motion is avail¬ 
able in the Supplemental Material). Apparently this pe¬ 
riodical motion of two bound vortices with opposite spins 
is induced by the polarization splitting in spinor conden¬ 
sates. It has been shown recently^ that similar spinor 
effects can evoke a uniform motion of bound polariton 
solitons in coherently driven microresonators. 


D. Scaling Laws for the Defect Density 

The number of vortices generated in the system outside 
of the MI region, during the mean-field evolution of an 
initially noisy, low-density state, is found to grow with 
the pump intensity, as shown in Fig. [8] The number of 
vortices is counted at a time after spatial coherence is 
established in the system, where the average polariton 
and reservoir density achieves a steady-state. Note that 
at very long-times there may be further recombination of 
vortex-antivortex pairs. 

By taking advantage of the universality of dynamics of 



FIG. 8: (Color online) (a) Same as in Fig. [T] under higher 
power excitation, (b) Variation of the number of vortices with 
pump power obtained numerically. The numerical results are 
averaged over ten different initial configurations at each value 
of pump intensity. The scaling with pump power is compared 
to the dependencies N& ~ P t h (1 — Pth/P) (green curve) and 
Nd ~ (P/Pth — 1) (magenta curve), see Eqs. (15) and |l6|. 


the system in the vicinity of the phase transition, we find 
approximate scaling laws governing the number of defects 
created during the transition, see solid line in Fig. [8](b) . 
The scaling laws in the case of exciton-polaritons have 
similar forms as t he on es obtained from the argument of 
Kibble and ZurekPI® 1 ! However, the dynamics of defect 
formation is different due to the fact that the initial state 
of the system is strongly out of equilibrium.^ While in 
the Kibble-Zurek mechanism the phase transition is as¬ 
sumed to begin in the initial state that is close to thermal 
equilibrium, here the white-noise initial state is domi¬ 
nated by fluctuations. In both cases, spontaneous sym¬ 
metry breaking occurs differently in separate regions of 
space which cannot communicate on how the symmetry 
is broken due to the finite timescale of the process. On 
the borders between these separate regions defects can 
appear in the form of domain walls, vortices, or more so¬ 
phisticated structures, depending on the dimensi onality 
of the system and the form of the order parameter ! 50 ! 51 ! 8 ^ 

In the case of dynamics described by 0 -<§ , start¬ 
ing from the initial low density white-noise state, the 
growth of polariton density initially occurs without es¬ 
tablishing any coherence (a movie showing the dynamics 
corresponding to the establishment of the state in Fig. [8ja 
can be found in the supplemental material). This corre¬ 
sponds to the first stage of the process, where there is 
practically no nonlinearity and no /c-dependence of the 
growth rate. The defects are created in the second stage, 
when fluctuations are suppressed due to nonlinear inter¬ 
actions. In this process, regions of ordered phase (or 
patches of a condensate) appear out of the initial dis¬ 
ordered strongly fluctuating phase. Due to spatial and 
polarization symmetry breaking, defects are created on 
the borders between these condensed regions. 

We assume that in the emerging ordered regions de¬ 
scribed above the relative phase between the two po¬ 
larization components is approximately equal to n and 
densities of |V’+| 2 and \i/j-\ 2 are similar. This assump¬ 
tion is in agreement with the numerical data presented 
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in Figs. [7]^c) and OH where the phase is equal to tt 
away from the eyelets, and is related to the fact that 
such configuration minimizes the energy in Eqs. (l|2|. 
We substitute i/j± = ±^ip to obtain 



dn R 

dt 


l ^ 2 V 2 oi + a 2 , , l2 \ , 

(--sr +!nn * + -2-1*1 r 

zfi 

+ y i rn R ~ 2r ) ^ - AxyV’, 

-{r R + \\if)n R + p. (10) 


when the condensate density is already near to its equi¬ 
librium value |V , o| 2 = (P — Pth) /r. Here we obtain 

N d ~ n 2 ~ 7o ~ P th (i - • (!5) 

Note that the above scaling does not have a power-law 
form, which is due to the fact that the transition is effec¬ 
tively nonlinear. 

In the opposite limit we assume that the vortices are 
formed when the condensate density is still very small 
I'i/’ol 2 = 0. In this case we obtain 


Next, we assume that the reservoir quickly adjusts to the 
change of 


UR r* + r#| 72 ' (11) 

Let us now assume that there exists a patch of approxi¬ 
mately constant polariton field -00 (or a condensate seed) 
at t = t\. We can describe small fluctuations around 7o 
by 

i/j = (ijj 0 + u(t)e ikx + v*(t)e ~ ikx ) e -vt+A(t- tl ) , (12) 


where /i and A are the chemical potential and the growth 
rate of the patch, and u, v represent the fluctuations. Ex¬ 
panding in Taylor series around ipo, we can rewrite 
(10b as 


ihij) 


ft 2 V 2 t 770 L_ 7| 2 \ 

2 m 2 y u. sa t ) 


T T /ia 




( 13 ) 

with the accuracy of the order of 0{\ip\ 4 ). Here 
9eff = (ai + « 2)(1 - ffirP/2r^)/2, 7o = (rP/T A )(l + 
r\ip 0 \ 2 /2T A ) - 2T, and n sat = 27 0 ry(r 2 P) where V A = 
Tfl + r |t/>o| 2 /2 is the effective exciton decay rate. The 
above equation has the same form as the one derived in 
Ref. [53], apart from the negligible dependence of g 0 g on 
P. The Bogoliubov-de Gennes modes are u, v* ~ e~ lul±t , 
and the mode frequencies 


~ = + M) 2 - (£fe + «/H (14) 

where £k = hk 2 /2m"fo, the saturation parameter a = 
l^o 1 2 Aw j an d P = ‘2gesF A /h r2 P. This spectrum has 
the property that modes with high momenta are strongly 
damped (the imaginary part of the frequency is nega¬ 
tive), in contrast to the initial linear dynamics where no 
fc-dependence of the imaginary part of the spectrum was 
present. We can define a characteristic momentum cutoff 
k for the modes that are strongly damped, which scales 
with the parameters as n ~ 7 0 ' . Any fluctuations with 
momenta higher than n will be suppressed, while fluctu¬ 
ations with lower momentum can form regions of ordered 
condensate phase. To estimate the scaling of the number 
of defects, two different limiting cases can be considered. 
In the first case, we assume that the vortices are formed 


N d ~ k 2 ~ 70 ~ H - 1. (16) 

-m th 

The two estimates [Eqs. and ( [T6| ] are compared to 
the numerical results in Fig. |sjb). The numerically ob¬ 
tained scaling appears to be intermediate between the 
two extreme cases. 


V. CONCLUSIONS 

In this paper we presented a comprehensive theoretical 
study of non-equilibrium dynamics of polariton conden¬ 
sates in incoherently pumped semiconductor microcav¬ 
ities. We have anticipated two different destabilization 
mechanisms that govern nonlinear dynamics of this sys¬ 
tem. The first arises when the polariton condensate has a 
strong feedback effect on the reservoir of incoherent “hot” 
polaritons. The second one is associated with the para¬ 
metric scattering in the presence of polarization splitting 
in a spinor condensate. Both mechanisms result in the 
formation of phase defects, i.e. vortices, triggered by the 
modulational instability of the homogeneous condensate. 

In the scalar case, we have shown that the presence 
of the incoherent reservoir can affect substantially both 
the vortex stability and their mutual collective dynamics. 
In particular this can lead to the formation of rotating 
dislocations or delocalized spiraling waves. 

In the spinor two-component condensate we have iden¬ 
tified the presence of topological defects, which take the 
form of half-vortex pairs, giving an “eyelet” structure in 
intensity and dipole type structure in the spin polariza¬ 
tion. 

In the case when the phase defects are formed in the 
dynamically (modulationally) stable regime, as a result 
of the condensate formation from an initial spatially and 
phase-incoherent state, we find that the defect density 
scales with the pumping rate in analogy to the Kibble- 
Zurek type scaling for non-equilibrium phase transitions. 
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Appendix A: Linear Stability Analysis 

The stability of the solutions can be checked using the 
standard approach of applying perturbations, in the form 
ip± | -t ( if)± + Sif >±), n R i->- n R + Sn R . Substitution 

into Eqs. 0;© and collecting terms linear in the small 
amplitudes dip- 1 and Sn R we have: 

= ci dip + + (A X y + a 2 ip*_ip + ) 6 i/j_ + ip+6il>+ 
+ a 2 V’-V’+ < ^/ ; - + (ffi + ihr) ip + Sn R , (Al) 
= c 2 ch/’- + (Aat + a^V’+V’-) dip+ + 

+ a 2 ip+ip-Sip *+ g 2 ip-5n R , (A2) 


= - (Tr + r\ip + \ 1 2 3 4 5 ) Sn R - 2 ma^te {ip+6il) + } , 

(A3) 

where we have defined: 

ci = gin R + ih(rn R -T) + 2a 1 \i/j + \ 2 +a 2 \ip-\ 2 , (A4) 
C 2 = g 2 n R - ihT + 2ai|'0_| 2 + a 2 \i/) + \ 2 . (A5) 


We consider perturbations in the polariton and reser¬ 
voir fields of the form 5ip ± = U±e i(.kx-wt) +v ^ e -i(kx-iu*t) 

and Sn R = w ( e Wx- u t) + e -i(kx-w*t)^ respectively. 
Here u± and v± are complex amplitudes, while w is a real 
amplitude, a; is a complex eigenvalue to be determined. 
Stability of the original solutions occurs if 3m {w} < 0 
such that the perturbation decays. 

Substitution of 6ip± and Sn R into Eqs. ( |A3[ ) and 
collection of terms oscillating as e~ lult and e~ lb2 1 yields 
a set of five coupled equations, which represent an eigen¬ 
value problem for u. In matrix form: 


/ c'r — Hu 
-aiipf 

A XY + a 2 il>*+4>- 
—a 2 ip*-ip+ 

\ - ihrn R ip* + 


OL\ll>% 

—d* — hui 
Ot2i>-1p+ 

-A XY + -a 2 ip + 4>*L 
—ihrn R ip+ 


A XY + a 2 il>lip + 
— a 2 ij)*_ ip+ 
c' 2 — hu> 

0 


a 2 il>-ip+ 

-Ayy - a 2 ip-il;* + 
aiip 2 . 

—c '2 — hu 
~ 0 


(31 + ihr)ip+ 

-( g 1 - ihr)ip + 
52^- 

-g 2 V~ 

-ih(T R + r\ip + \ 2 ) - hw 


\ 

/ u+ \ 
V+ 


U- 

/ 

V- 

\ W / 


= 0 , (A6) 


where = Ci + and c ' 2 = c 2 + AT— hg. 

The imaginary and real parts of the perturbation spec¬ 
tra, E = hu, are shown in Fig. |9](a,b). The energies of 
the stationary states hfj, are illustrated by thin horizontal 


lines in Fig. [9](a). Where the imaginary parts are posi¬ 
tive, thick curves denote modulational instability at the 
given wavevectors k. The stationary state labelled U 2 is 
unstable to fluctuations with k = 0. 
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